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1 Abstract 

We present the results of Suzaku observation of the radio halo cluster Abell 2319. The metal abundance 
^ . in the central cool region is found to be higher than the surrounding region, which was not resolved in the 

former studies. We confirm that the line-of-sight velocities of the intracluster medium in the observed region 
are consistent with those of the member galaxies of entire A2319 and A2319A subgroup for the first time, 
though any velocity difference within the region is not detected. On the other hand, we do not find any 
signs of gas motion relevant to A2319B subgroup. Hard X-ray emission from the cluster is clearly detected, 
but its spectrum is likely thermal. Assuming a simple single temperature model for the thermal component, 
we find that the upper limit of the non-thermal inverse Compton component becomes 2.6 x 10 -11 erg s _1 
. cm -2 in the 10-40 keV band, which means that the lower limit of the magnetic field is 0.19 fiG with the 

Qh radio spectral index 0.92. Although the results slightly depend on the detailed spectral modeling, it is 

q , robust that the upper limit of the power-law component flux and lower limit of the magnetic field strength 

become ~3x 10 -11 erg s _1 cm -2 and ~ 0.2/iG, respectively. Considering the lack of a significant amount of 
very hot (~ 20 keV) gas and the strong bulk flow motion, it is more likely that the relativistic non-thermal 
electrons responsible for the radio halo are accelerated through the intracluster turbulence rather than the 
shocks. 
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1/-} i 1. Introduction processes. However, it is not so easy to determine clus- 

■ ter magnetic field strength and structures in an obser- 

Diffuse non-thermal synchrotron radio emission is found vational way. Faraday rotation measure observations of 
. in a significant fraction of galaxy clusters (Giovannini & polarized radio sources in and/or behind clusters are of- 
' Feretti 2000; Kempner & Sarazin 2001), which indicates ten used for this purpose (Clarke et al. 2001; Vogt & 
ON . that there exist both the relativistic electrons and mag- Ensslin 2003; Govoni et al. 2006). However, the results 
netic field as well as the thermal intracluster medium strongly depend on the magnetic field structures them- 
(ICM) in the intracluster space. Some of such diffuse selves. Moreover, it can be applied for very limited re- 
radio sources are located near the center and cover the gions where we have suitable polarized radio sources by 
cluster entirely, which are called "radio halos". Others chance, though this will be possibly resolved by using 
are located at periphery and called "radio relics" . cosmic microwave background (CMB) as the polarized 

Although the origin of these non-thermal electrons is sources (Ohno et al. 2003; Takizawa 2008). Therefore, it 
still unclear, some connections of radio halos and relics is very important to estimate the magnetic field strength 
with dynamical motion of ICM are reported. For example, through another independent method, 
they are rarely found in cool core clusters (Giovannini et GeV electrons in radio halos and relics are expected to 
al. 1999). Radio luminosity has strong correlation with emit non-thermal hard X-ray via inverse Compton pro- 
X-ray one and ICM temperature (Liang et al. 2000). In cess of CMB photons. Comparing the synchrotron radio 
addition, Buote (2001) shows the correlation of the radio flux and inverse Compton hard X-ray one (or its upper 
power with asymmetry of the X-ray image. limit), we are able to estimate volume- averaged magnetic 

It is likely that such non-thermal electrons are ac- field strength (or its lower limit). However, it is very 
celerated via shocks (Sarazin 1999; Takizawa & Naito difficult to detect such components because of a large 
2000; Totani & Kitayama 2000; Miniati et al. 2001; Ryu amount of thermal emission from the ICM. Although a 
et al. 2003; Takizawa 2003; Inoue et al. 2005) and/or lot of efforts have been made to detect the non-thermal 
turbulence (Roland 1981; Schlickeiser et al. 1987; Blasi inverse Compton component from clusters of galaxies, 
2000; Ohno et al. 2002; Fujita et al. 2003; Brunetti et the situation is still unclear. For example, detection of 
al. 2004) in the ICM associated with cluster mergers, non-thermal hard X-ray from the Coma cluster were re- 
Magnetic field plays a crucial role in these acceleration ported from the Beppo-SAX PDS (Fusco-Femiano et al. 



2 

1999; Fusco-Femiano et al. 2005) and RXTE (Rephaeli et 
al. 1999; Rephaeli&Gruber 2002) data, but their reliabil- 
ity is still controversial (Rossetti, & Molendi 2004; Fusco- 
Femiano et al. 2007). Recently, in addition, Wik et al. 
(2009) show the lower limit that conflicts with the former 
detection reports through combined analysis of Suzaku 
and XMM. Therefore, firm detection or reliable upper lim- 
its of non-thermal components in the hard X-ray regime 
are highly desired by independent instruments. 

Abell 2319 is one of the most well-known examples of 
merging clusters with a giant radio halo. Two subgroups, 
A2319A and A2319B, are recognized in radial velocity dis- 
tribution of the member galaxies, which suggests that the 
merger axis is nearly along the line-of-sight and that the 
velocity difference between them is almost 3000 km s _1 
(Oegerle et al. 1995). Markevitch (1996) studies tem- 
perature structures in a relatively large scale and finds 
that the temperature in the region corresponding to the 
A2319B is lower. Chandra observations reveal inhomo- 
geneous temperature structures and a cold front in the 
central region, and show that the position of the X-ray 
center is different than that of the cD galaxy (O'Hara et 
al. 2004; Govoni et al. 2004). Thus, it is obvious that the 
cluster is not dynamically relaxed. A2319 has a giant ra- 
dio halo, and the radio and thermal X-ray distributions 
are quite similar to each other (Feretti et al. 1997; Govoni 
et al. 2001). Observations in the hard X-ray band are 
performed by Beppo-SAX (Molendi et al. 1999), RXTE 
(Gruber&Rephaeli 2002), and Swift-BAT (Ajello et al. 
2009), any of which do not report firm detection of non- 
thermal components. 

The hard X-ray detector (HXD) on board Suzaku 
(Mitsuda et al. 2007) is superior in order to investigate 
hard X-ray properties of galaxy clusters because of its low 
detector background and narrow field of view (Kokubun 
et al. 2007; Takahashi et al. 2007). Indeed, improved con- 
straint for non-thermal components has been obtained 
for several clusters (Fujita et al. 2008; Kawano et al. 
2009; Nakazawa et al. 2009; Wik et al. 2009), and very 
hot components (~ 20 keV) were found in RXJ1347.5- 
1145 (Ota et al. 2008) and A3667 (Nakazawa et al. 2009). 
In addition, high spectral resolution of the XIS (Koyama 
et al. 2007) on board Suzaku makes it possible to constrain 
the bulk flow motion of the ICM that is likely relevant to 
the particle acceleration processes (Ota et al. 2007; Fujita 
et al. 2008). 

In this paper, we present Suzaku observation of the 
Abell 2319 cluster to investigate dynamical status of the 
ICM and hard X-ray properties, which enables us to un- 
derstand particle acceleration processes and the origin of 
the radio halo. The canonical cosmological parameters 
of H = 70 Mpc -1 km s" 1 , ^ = 0.3 and A = 0.7 are 
used in this paper. Because the cluster mean redshift is 
z = 0.0557, 1' corresponds to 62 kpc. Unless otherwise 
stated, all uncertainties are given at the 90% confidence 
level. 

The rest of this paper is organized as follows. In §2 
we describe the observation and data reduction. In §3 we 
present the spectral analysis results. In §4 we discuss the 
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Fig. 1. A ROSAT PSPC image of the Abell 2319 overlaid 
with field of view of Suzaku XIS CCDs (blue), and that of 
HXD PIN (magenta) in which the throughput of the fme-col- 
limator becomes 50%. An approximate position of the sub- 
group A2319B is also represented by dotted yellow circle. 

results and their implications. In §5 we summarize the 
results. 

2. Observation and Data Reduction 

We observed the central region of Abell 2319 with 
Suzaku on 2006 October 27-30 for an exposure time of 100 
ks. The field of view (FOV) of Suzaku XIS, and that of 
HXD PIN in which the throughput of the fine-collimator 
becomes 50 % are shown in a ROSAT PSPC image (fig- 
ure 1). An approximate position of the A2319B subgroup 
is also shown. The observation was performed at HXD 
nominal pointing. As a result, most of the A2319B is not 
within the XIS FOV. The XIS was operated in the normal 
full- frame clocking mode. The edit mode was 3x3 and 
5x5, and we used the combined data of both modes. 
A half of the 64 PINs of the HXD were operated in the 
nominal bias voltage of 500 V, and the others were done in 
that of 400 V. All the data was processed with the Suzaku 
pipeline processing of version 2.0.6.13. We employed the 
calibration data files of 20081009. 

The XIS data were processed through standard crite- 
ria as follows. Events with a GRADE of 0,2,3,4,6 and 
STATUS with 0:524287 were extracted. We exclude the 
data obtained at the South Atlantic Anomaly (SAA), 
within 436 sec after the passage of SAA, and at the low 
elevation angles from the Earth rim of < 5° and the sun-lit 
Earth rim of < 20°. As a result, effective exposure time 
becomes 99.5 ks for XIS. Non X-ray background (NXB) 
spectra and images of XIS were generated using the ftool 
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Fig. 2. XIS image of the A2319 cluster in 0.5-8.0 keV band 
combined from the FI CCD images. The images were cor- 
rected for exposure and vignetting effects after subtracting 
NXB, and smoothed by a Gaussian kernel with a = 0.3 ar- 
cmin. 

"xisnxbgen" (Tawa et al. 2008). Figure 2 represents the 
0.5-8.0 keV XIS image combined from those of the front 
illuminated (FI) CCDs (XIS0, XIS2, XIS3). The images 
were corrected for exposure and vignetting effects after 
subtracting NXB, and smoothed by a Gaussian kernel 
with a = 0.3 arcmin. 

The HXD data were also processed in a standard way. 
We exclude the data obtained at SAA, within 500 sec after 
the passage of SAA, within 180 sec before entering SAA, 
at the low elevation angles from the Earth rim of < 5°, 
and at the location where the Cut- Off- Rigidity (COR) is 
lower than 6 GV. The resultant effective exposure time is 
94.0 ks. As the NXB of the PIN, we used a public NXB 
model provided by the HXD team in the Suzaku website. 
The version of the model is "METHOD=LCFITDT" , or 
"tuned" (Fukazawa et al. 2009). The CXB level is esti- 
mated in the same way as Nakazawa et al. (2009) from 
the Lockman hole observation (Suzaku observation ID, 
101002010). The detailed procedure is described in ap- 
pendix 2 of Nakazawa et al. (2009). We defined the photon 
flux model as N(E) = 8.7 x 10~ 4 x E' 1 ' 29 x exp(-£/40.0) 
in photons cm -2 s _1 keV -1 FOV -1 , where E is the pho- 
ton energy in keV. 

3. Spectral Analysis 

For spectral analysis of the XIS data, rmf and arf files 
are generated using the ftool "xisrmfgen" and "xissimarf- 
gen" , respectively. The ROSAT PSPC image of the A2319 
is used to make arf files, considering that the XIS FOV 
does not cover the cluster ICM emission entirely and that 
Suzaku has only moderate spatial resolution. The clus- 



ter ICM emission in the XIS FOV is so bright that the 
background components from astrophysical origin such as 
CXB and Galactic diffuse components are safely negligi- 
ble. 

We use an rmf file provided by the HXD team in the 
CALDB for PIN spectral analysis. A PIN arf file are made 
in the same way of Nakazawa et al. (2009); we calculate 
it by convolving the point source arf assuming that the 
spatial distribution of the emission is the same as that of 
the ROSAT PSPC image. Detailed description is written 
in subsection 2.4 of Nakazawa et al. (2009). 

3.1. Temperature and Abundance Structures 

Temperature and abundance structures can provide us 
with basic and crucial information about the ICM dynam- 
ics and its origin. To explore this, we divide the region 
observed with XIS into 15 regions presented in figure 3 and 
perform spectral analysis for each region. The sizes of the 
divided region are 3' x 3' and 6' x 6' for the cluster cen- 
ter and surrounding regions, respectively. FI CCD spectra 
and related response files for each region are summed, and 
FI and back illuminated (BI) CCD spectra are fitted si- 
multaneously. For the spectral fitting, we use the energy 
band of 0.6 - 10.0 keV and 0.6 - 8.0 keV for XIS FI and 
BI, respectively, though the energy band near the Si edge 
(1.7-1.9 keV) is ignored to avoid uncalibrated structures. 
Each spectrum is fitted by the photoabsorbed single tem- 
perature APEC model (WABS x APEC), assuming that 
N H = 7.93 x 10 20 cm" 2 (Dickey & Lockman 1990) and 
that redshift=0.0557 (Struble & Rood 1999). The fitting 
results are presented in table 1. In general, each data 
are fitted well with the above-mentioned model. Figure 
4 and 5 show the resultant temperature and metal abun- 
dance for the regions presented in figure 3, respectively. 
The central south-east regions (region 2, 3, and 4) tend 
to have lower temperature and higher abundance. 

In order to confirm this trend of temperature and metal 
abundance distribution, we determined the X-ray spec- 
tra in annuli whose center is located at the X-ray peak 
(290.3000, +43.9458). The details of the spectral fitting 
procedure are similar in 2D analysis presented above. 
Figure 6 and 7 show the radial (not deprojected) profile of 
the temperature and abundance, respectively. It is clear 
that the ICM in the central region has lower temperature 
and higher abundance. 

3.2. Line- of- sight Velocities of the ICM 

Line-of-sight velocities of the ICM contain more direct 
information about the cluster dynamics than X-ray sur- 
face brightness and temperature structures. In general, 
however, more photons are necessary to determine the 
center of the Doppler-shifted lines with meaningful ac- 
curacy through the spectral fit of the XIS data than to 
determine the temperature and metal abundance, which 
means that we have to use larger regions. Thus, we divide 
XIS FOV into 11 regions shown in figure 8 and perform 
spectral analysis for each region. The sizes of a region are 
3' x 3' and 6' x 9' for the central regions (from 1 to 6) 
and the outer parts (from 7 to 11), respectively. In order 
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Table 1. Best fit parameters for the XIS spectra of the re- 
gions presented in figure 3. 



region 


kT 


Z 






(keV) 


(Ze) 


1 


ii -°-0.4 


0.27^;- 


977.3/892 


2 


Q 7+ - 2 

y -' -0.2 


o si +0 - 02 

u - o± -0.02 


1540.0/1389 


3 


Q fi+ - 2 

y -°-0.2 


n 9 o+0.03 


967.9/982 


4 


9.6t° ; 2 2 


u - zy -0.03 


1198.3/1132 


5 


I0.li° ; 2 2 


U.O4_ q 2 


1775.7/1610 


6 


I0.5i° ; 2 2 


24 +0 - 03 

u - z ^-0.03 


1117.3/990 


7 


11 7 +0 - 4 


n 9O+0.03 
u - zc) -0.03 


1162.3/1079 


8 


ID 1+ - 3 


90+O.O3 
u - zo -0.03 


953.6/902 


9 




u - zo -0.03 


1219.0/1127 


10 


io.it° : 2 2 


24+ - 03 


1012.3/934 


11 


10.4lg;i 


n 90+O.O5 
u - zo -0.05 


501.3/449 


12 


io.4i8;i 




420.2/445 


13 


8.7^1 


20 +0 - 06 

U ' ZU -0.06 


283.0/282 


14 


8.9t° : 3 3 


21 +0 - 04 


623.1/551 


15 


9.7t° : 6 6 


21 +0 - 07 
U ' Zi -0.07 


362.8/251 
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Fig. 4. Temperature for the regions presented in figure 3. 
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Fig. 5. Same as figure 4, but for the metal abundance. 



Fig. 3. Region numbers used in two-dimensional tempera- 
ture and abundance distribution measurement in figure 4 and 
5. 
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Fig. 6. Projected radial temperature profile. 
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Fig. 7. Same as figure 6, but for the metal abundance. 

to determine ICM line-of-sight velocities, we mainly use 
Doppler-shifted He like Fe Ka lines (6.679 keV) and H 
like Fe Ka line (6.964 keV). Therefore, if the ICM motion 
is ~ 1000 km s _1 , we should resolve the energy shift of 
only ~ 22 eV. This is a challenging task considering that 
energy resolution of XIS is 130 eV (FWHM). 

To measure an energy scale of XIS accurately, we make 
a gain correction in a similar way of Fujita et al. (2008) 
with Mn Ka lines of the calibration sources on each XIS 
sensor. We made spectra of the calibration source region 
for each XIS sensor, and fit them with APEC plus two 
Gaussians assuming that the abundance and redshift of 
APEC are equal zero. A gain correction factor is defined 
as 



gam 



£ fit (Mn-Ka) 
E (Mn-Ka) ' 



(1) 



where ^t(Mn-Ka) and £?o(Mn-Ka) are the energy value 
obtained from the fitting results and expected one (5.895 
keV), respectively. We evaluate a gain correction factor 
/gain for each XIS sensor. Typically, we find | /gain — 1| 
to be ~ 0.0005, which corresponds to be ~ 150 km s _1 . 
Then, the corrected redshift becomes 



^cor — /gain(^fit H~ 1) 1? 



(2) 

where z cor and Zf& are the corrected redshift and that 
obtained from the fitting results, respectively. 

Although the above-mentioned correction is obviously 
valid for the regions near the calibration sources, the gain 
could depend upon the position on the same CCD chip ow- 
ing to the charge transfer inefficiency (CTI). As for this, 
an accuracy of 0.2% was reported over the CCD chips be- 
fore spaced-row charge injection (SCI) was adopted (Ota 
et al. 2007). Matsumoto et al. (2007) 1 reported that CTI 
was almost zero just after the SCI was adopted. Since 
our observation was held at the timing very near to this 
report, we conclude that the systematic errors of the en- 
ergy scale because of CTI are estimated to be 0.2%, which 
corresponds to be ~ 600 km s _1 . 




Fig. 8. Region numbers used in measurement 
line-of-sight velocity of the ICM in figure 9. 



of the 



We fit spectra of each XIS sensor with APEC model in 
the energy band of 5.0 - 10.0 keV. The resultant values 
of Zfit are corrected into z cor in a way of equation (2), 
which are converted into line-of-sight velocities. Then, 
we calculate arithmetic means of the results of all four 
XIS sensors. Figure 9 shows line-of-sight velocities of the 
ICM for the regions presented in figure 8, where red , 
green, and blue solid lines represent the mean line-of-sight 
velocities of the member galaxies for the entire A2319, 
A2319A subgroup, and A2319B subgroup, respectively. 
The error bars stand for only statistical ones. Basically, 
the obtained ICM velocities lie between those of entire 
A2319 and A2319A. No significant velocity difference is 
detected within the observed region. Taking account of 
the systematic errors because of CTI, we confirm that the 
observed ICM velocities are consistent with those of entire 
A2319 and A2319A subgroup, and do not find any signs 
of ICM motion associated with A2319B subgroup in the 
observed region. 

As mentioned before, a cold front is found in the south- 
east of the X-ray center by Chandra observations (O'Hara 
et al. 2004; Govoni et al. 2004). In addition, we find that 
the region in the north of the cold front has lower temper- 
ature and higher abundance. This may mean that there 
exits a cold core originated from the subcluster that in- 
falls in the past. Therefore, it is possible that ICM in this 
region has negligible relative velocity to the surrounding 
ICM. In order to examine this, we analyze spectra of the 
regions presented in figure 10 and measure the line-of-sight 
velocities in a similar way. The results are shown in figure 
11. Again, no significant velocity difference is detected. 
The obtained ICM velocities lie between those of entire 
A2319 and A2319A. It is also clear that there are no signs 
of ICM motion related with A2319B subgroup. The size 
of Region 1 in figure 10 is approximately 3' x 5.5'. Taking 
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Fig. 9. Line-of-sight velocities of the ICM for the regions pre- 
sented in figure 8. Red , green, and blue solid lines show the 
mean line-of-sight velocities of the member galaxies for the 
entire A2319, A2319A subgroup, and A2319B subgroup, re- 
spectively. The error bars stand for only statistical ones. 



Fig. 11. Same as figure 9, but for the regions presented in 
figure 10. 




Fig. 10. Regions used in measurement of the line-of-sight ve- 
locity around the cold front in figure 11. 

account of the X-ray telescope response, we estimate that 
~ 30% of photons detected in the Region 1 are originated 
from Region 2, and less than 10%, vice versa, respectively. 
Therefore, this contamination does not affect our results 
significantly. 



3.3. Hard X-ray Properties 
thermal Components 



and Constraint of Non- 



The PIN spectrum of this observation is shown in Figure 
12, where black, red, and green crosses represent the data, 
NXB model, and residual signals (data-NXB), respec- 
tively. A typical CXB model spectrum is also shown in 
blue crosses. Hard X-ray emission from A2319 is clearly 
detected in the energy below ~ 40 keV. Therefore, we use 




Energy (keV) 

Fig. 12. The PIN spectrum of this observation. Black, red, 
and green crosses show the spectrum of the data, NXB model, 
and residual signals (data-NXB), respectively. A typical CXB 
model spectrum is also presented in blue crosses. 

the PIN spectrum in the energy band of 13.0-40.0 keV 
hereafter. We also check the GSO spectrum and confirm 
that the data are consistent with NXB. Thus, we do not 
use them. 

First, we fit the PIN spectrum alone with the APEC 
or power-law model in the energy range of 13.0-40.0 keV. 
The metal abundance and redshift in the APEC model 
are fixed to be 0.3 and 0.0557, respectively, because it 
is very difficult to determine these parameters with only 
PIN data. The results of the fit are presented in table 2. 
It is obvious that the spectrum is better represented by a 
thermal APEC model than a non-thermal power-law one. 
The temperature determined from the PIN spectrum is 
similar to the values in figure 4 determined from the XIS 
data. 

Second, we perform joint spectral analysis of PIN and 
XIS to investigate hard X-ray properties. It should be 
noted that soft X-ray spectral information is essential to 
constrain spectral components characteristic in the hard 
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Table 2. Fitting results of the PIN spectrum 



model 


T or fcT(keV) 




PL 


3 l+°- L 
°- i -o.i 


81.7/69 


lkT 


in q +0 - 9 


67.7/69 



X-ray band such as non-thermal power-law and/or very 
hot thermal ones. We use the XIS spectrum in the en- 
ergy band of 2.0-10.0 keV and 2.0-8.0 keV for FI and BI 
CCDs, respectively. FI CCD spectra and associated re- 
sponse files are summed, and FI, BI, and PIN spectra are 
fitted simultaneously. Systematic errors between XIS and 
PIN normalization are taken into account properly. We 
tried to fit XIS+PIN spectra in a similar way but for using 
XIS spectra also below 2 keV. The fitting results are not 
acceptable if we adopt common Nh for FI and BI CCDs, 
because significant residuals inconsistent between FI and 
BI CCDs are found in the soft band. Then, we tried to 
fit the spectra for each XIS sensor and PIN simultane- 
ously allowing that each XIS has a different Nh value. 
This improves the fitting results significantly. We suspect 
that this is possibly because the contamination correction 
of arf files are not very good in case of regions as large 
as almost entire XIS FOV, though we have no definitive 
idea about this problem at present. In addition, it is pos- 
sible that this kind of systematic errors become appar- 
ent because the spectra have sufficiently good statistics. 
The difference of Nh among four XIS sensors is typically 
~ 1 — 3 x 10 20 cm -2 , which is so small that it does not se- 
riously affect the results of the spectral fit above 2.0 keV. 
Thus, we decide not to use XIS spectra below 2.0 keV to 
avoid this unsolved issue. 

Since we are interested in both non-thermal inverse 
Compton and very hot components, spectral models used 
here are a single temperature APEC model (lfcT), single 
temperature APEC plus power-law model (lfcT + PL), 
two temperature APEC model (2fcT), and two temper- 
ature APEC plus power-law model (2kT + PL). Again, 
we fix both N H = 7.93 x 10 20 cm -2 and redshift=0.0557 
as in §3.1. In the 2kT and 2kT + PL models, a relative 
normalization ratio between two APEC models is allowed 
not to be common between XIS and PIN, to compensate 
small (~ 10%) uncertainty in the arf generation. It is 
very difficult to constrain power-law components in the 
XIS energy range because of dominant presence of thermal 
emission, which could degrade the HXD fitting results. In 
order to avoid this, thus, a relative normalization ratio 
between APEC and power-law models is also allowed not 
to be common between XIS and PIN. A photon index 
of the power-law component is also fixed assuming that 
it is emitted via inverse Compton process of CMB from 
the same electron population that radiates synchrotron 
radio. Considering that typical magnetic field strength of 
the intracluster space is ~ /xG, electrons that radiates the 
inverse Compton hard X-ray of 20-40 keV emit relatively 
low frequency synchrotron radio. We adopt 1.92 or 2.4 
as the photon index of the power-law component because 
0.92 (408-610MHz) and 1.4 (26-610MHz) are reported as 



Table 3. Best fit parameters of the lkT and 2kT models for 
XIS and PIN spectra. 



lkT 2kT 



fcT low (keV) 


9 7+ - 1 

y - ' -0.1 


6-5^5 


fcT high (keV) 






Z(Z Q ) 


26+ - 01 
u - zo -o.oi 


0.28±g;S 


Mow, XIS 


o.i7lS:82 


8.5 + i'i x 10" 2 


Mow, PIN 


.i8ir i 


14+' - 04 

U ' i4 -0.04 


A^hi g h,xis 




9.2+^ x l(r 2 


Miigh,PIN 




6A±\ A 3 x 10" 2 


X 2 /d.o.f. 


3215.6/2978 


3136.6/2975 



Table 4. Same as table 3, but for the lkT + PL models, 
where a photon index of the power low component is fixed to 
be 1.92 or 2.4. 



lkT + PL(1.92) lkT + PL(2.4) 



A:T low (keV) 


9 7+ - 1 

y - ' -0.1 


9 7+ - 1 

y - ' -0.1 


Z(Z Q ) 


26+ - 01 


26+ - 01 
u - ZD -o.oi 


Mow, XIS 


17+ - 00 

u ' i ' -0.00 


17+ - 00 
u - i ' -0.00 


Mow, PIN 


16+ - 02 

U - 1D -0.01 


u - lo -0.02 


r PL 


1.92(fixed) 


2.4(fixed) 


MpL,XIS 


0.0±g;£ x 10- 3 


O.Ot^ x 10- 3 


MpL,PIN 


3.8t 2 | x 10- 3 


2.2ti s 7 x 10~ 2 


X 2 /d.o.f. 


3210.9/2976 


3211.2/2976 



radio spectral index (Erickson et al. 1978; Harris & Miley 
1978; Feretti et al. 1997). 

Figure 13, 14, and 15 show the spectral fit results for 
the lkT, 2kT, and lkT + PL model with the photon in- 
dex 1.92, respectively. Best fit parameters of the lkT and 
2kT models are listed in table 3. Although the reduced 
chi square is slightly improved in the 2kT model, the lkT 
model gives us sufficiently acceptable results. Thus, it 
is possible that there exits a little very hot component 
(~ 16keV), which does not have to be needed. The de- 
tailed fitting results of the lkT+PL and 2/cT+PL models 
are also presented in table 4 and 5, respectively. Clearly, 
adding a power-law component does not improve the fit- 
ting results significantly in both the lkT and 2kT models. 

Our spectral modeling for thermal components are 
somewhat too simple considering that Chandra observa- 
tions (O'Hara et al. 2004; Govoni et al. 2004) show multi 
temperature structures with a certain temperature range. 
In general, if multi-temperature spectra are forced to be 
fitted with a simple lkT ^ PL model, the power-law com- 
ponent tends to be overestimated. Thus, the derived up- 
per limit of the power-law component can be regard as 
relatively conservative one. In case of the 2kT modeling, 
on the other hand, the situation could be more compli- 
cated. The higher energy part of the spectra can be rep- 
resented by both higher temperature and/or power-law 
components. It is not trivial how this increase of the de- 
gree of the freedom of the model affects the actual fitting 
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Table 5. Same as table 4, but for the 2kT + PL models. 



2kT + PL(1.92) 2kT + PL(2A) 



rvJ - low V v / 


6 5 + "-2 

u * u — 0.6 




6 5+^ 

u-u — 0.6 




kTh\<rh(keV) 

,v -*- nigii y- 1 *-^ » ) 


15.7+?-? 

u 1 — 1.1 




15.7+?-? 

- L ^* ' — 1.1 




V W / 


0.28+S-2J 

—0.01 




0.28+°°j 

—0.01 




^V"low,XIS 


8.6 + j-i x 10- 


-2 


8.6 + {-j x 10- 


-2 


Mow, PIN 


U ' iO -0.04 




u - iO -0.04 




^high,XIS 


9.2±^| x 10" 


-2 


9.2t^ x 10" 


-2 


Miigh,PIN 


6.3^ x 10" 


-2 


6.3±£ : | x 10" 


-2 


r PL 


1.92 (fixed) 




2.4(fixed) 




^VpL,XIS 


0.0±g;^ x 10" 


3 


0.0±g;g x 10" 


-4 


^V"PL,PIN 


0.0j:J : g x 10" 


-2 


0.0±g : g x 10" 


-2 


x7d-o-f- 


3136.7/2973 




3136.7/2973 






Energy (keV) 

Fig. 13. The wide band spectrum of the A2319 fitted with 
the lkT model. Black, red, and green crosses show the spec- 
trum of XIS-FI, XIS-BI, and HXD-PIN, respectively. The 
best fit model is presented in dashed histograms. 




5 10 20 

Energy (keV) 



Fig. 14. Same as figure 13, but for the 2kT model. Two 
thermal components are independently presented as dashed 
histograms. 




5 10 20 

Energy (keV) 



Fig. 15. Same as figure 13, but for the lkT-\-PL model with 
a photon index 1.92. Thermal and power-law components are 
independently presented as dashed histograms for HXD. Only 
thermal component is presented for XIS. 

results. 

In order to constrain an upper limit of a non-thermal 
power-law component, not only statistical errors but also 
systematic ones of both CXB and NXB have to be taken 
into account properly. It is well-known that the CXB 
fluctuations can be modeled as <tcxb/^cxb cx Sl" ' 5 ^ 25 , 
where Vt e and S c are the effective solid angle and up- 
per cutoff flux of a point source, respectively. From the 
HEAO-1 A2 results, ct C xb//cxb = 2.8% with n e = 15.8 
deg 2 and S c = 8 x lO" 11 erg s" 1 cm" 2 (Shafer 1983). 
Considering a rough estimate of NXB uncertainties of PIN 
(~5xl0 -12 erg s _1 cm -2 ), we adopt conservative upper 
cutoff flux S c ~8x 10- 12 erg s" 1 cm" 2 . Thus, the CXB 
fluctuation in the HXD PIN FOV is expected to be 18% 
at the 90% confidence level because of Q e = 0.32 deg 2 . 
Fukazawa et al. (2009) reports that the reproducibility of 
blank sky observations separated into 10 ks exposures give 
distribution of 5.7% at the 90% confidence level, where 
contribution from the statistical error of typically 3.3 % or 
larger and the effect of CXB fluctuation (1.3% of the total 
background) are included. This means that the NXB sys- 
tematic uncertainty is 4.5% (= - (3.3) 2 - (1.3) 2 ) 
at the 90% confidence level. 

In the joint spectrum analysis of PIN and XIS, NXB and 
CXB components are fluctuated at the 90 % confidence 
level of the systematic uncertainty mentioned- above (4.5% 
and 18% for NXB and CXB, respectively) . This causes 
the changes of the best-fit parameters in the fits and gives 
us the systematic errors. Including both the statistical 
errors and systematic ones of CXB and NXB, we derive 
an upper limit of a power-law component in 10-40 keV at 
the 90 % confidence level for each model. The results are 
presented in table 6, where the former and latter errors 
of flux are statistical and systematic ones, respectively. 
Although the results slightly depend on spectral modeling, 
all models give us the upper limit of ~ 3 x 10 -11 erg s _1 
cm -2 . 
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Table 6. Flux of a power-law component in 10-40 keV with 
statistical and systematic errors and its upper limit at the 90 
% confidence level for each model. 



model 



flux 

(erg s _1 cm -2 ) 



upper limit 

(erg s~ 



cm 



1/cTH 


-PL(1.92) 


-1 -1+0.8+1.3 
0. 8-1.1 


X 10" 


ii 


< 2.6 x 10- 


ii 


lkT-\ 


-PL(2A) 


-1 r + 1.2+1.9 

J-' -!. 2-2.1 


x 10- 


11 


<3.8x 10- 


n 


2/cTH 


-PL(im) 


n n +3.7+2.2 

u - u -o. 0-0.0 


x 10- 


11 


<4.3x 10- 


n 


2/cTH 


-PL(2A) 


n n +2.2+3.9 
u ' u -0. 0-0.0 


x lO- 


11 


<4.5 x lO- 


li 



4. Discussion 

4-1. Temperature and Abundance Structures 

Complex temperature structures are obtained in the 
central region of the Abell 2319 by Chandra observations 
(O'Hara et al. 2004; Govoni et al. 2004). Basically, our 
results are consistent with theirs though FOV of both ob- 
servations are not perfectly overlapped with each other. 
There are high temperature regions in the north east of 
the X-ray center (region 1 and 7 in figure 3), and regions 
in the north-west of the cold front (region 2, 3, and 4 
in figure 3) have lower temperature. However, very high 
temperature (~ 15 keV) components in the Chandra data 
are not found in our 2D spectral analysis of the XIS data 
alone. On the other hand, the 2kT model fitted to the 
XIS+PIN spectrum shows a strong 15 keV component, 
but this is likely an artifact from introducing a strong 6.5 
keV component for a cluster with average temperature of 
9.7 keV. Since the spectral fit with only HXD prefers a 
temperature less than 11.8 keV, and the XIS+PIN wide 
band spectrum is well fitted by a 9.7 keV single temper- 
ature thermal model, we conclude that there is no strong 
evidence of a bright very hot component, as high as 15 
keV in the Suzaku data. It should be noted that tem- 
perature determination of very hot gas with temperature 
more than 10 keV is difficult for Chandra CCD without 
any sensitivity above 10 keV. It is also possible that our 2D 
XIS results are significantly spatially smoothed because of 
the moderate spatial resolution of Suzaku. 

Markevitch (1996) studied a larger scale temperature 
structure of the A2319 with ASCA. It reveals that the 
A2319b region and outer parts of the cluster have lower 
temperature, which are without FOV of XIS but partly 
within that of PIN. Our wide band spectral analysis of 
the 2kT model shows that a relative ratio of normaliza- 
tion between two APEC models in XIS and PIN are signif- 
icantly different than each other. The results tell us that 
the lower temperature gas is more dominant in the PIN 
spectrum than XIS, which is not surprising and naturally 
expected from the above-mentioned ASCA results. 

Our 2D and projected radial spectrum analysis show 
that metal abundance in the central cool region is clearly 
higher than the surrounding region (see figure 5 and 7). 
Molendi et al. (1999) performed similar analysis with 
Beppo-SAX, and shows that the metal abundance distri- 
bution is consistent with homogeneous one within statis- 



tical errors. Because of large effective area and high sen- 
sitivity of XIS, metal abundance inhomogeneity in A2319 
is clearly detected for the first time with Suzaku. 

4.2. Bulk Flow Motion in the ICM 

As written before, though no significant line-of-sight ve- 
locity difference within the observed region is detected, we 
find that their values are consistent with those of entire 
A2319 and A2319A subgroup for the first time. In figure 
9, the largest velocity difference is found between region 
4 and 8, which provides us with the value 940 1 ^31 km 
s _1 . The sound velocity of 10 keV ICM is - 1700 km s _1 . 
Thus, it is probable that the internal motion of the ICM in 
the observed region is subsonic, though possibly it is close 
to the sound velocity. On the other hand, the difference 
between the observed region and A2319b is nearly 3000 
km s _1 . This means that super sonic ICM collision is ex- 
pected if the ICM associated with the A2319b has the ve- 
locity similar with the A2319b member galaxies. Another 
possibility is that the A2319b group is not gravitationally 
bound to the A2319 but located along the line-of-sight by 
chance. Please note that this possibility is not negligible 
from orbital motion analysis with a two-body model in 
Oegerle et al. (1995). 

4.3. Constraint of the Nonthermal Emission and 
Magnetic Field Strength 

Let us constrain the magnetic field strength from the 
combined analysis of radio synchrotron flux and an inverse 
Compton hard X-ray upper limit (Rybicki & Lightman 
1979). Typical energy hv' of photons scattered via in- 
verse Compton processes by electrons with energy jm e c 2 
is hv' ~47 2 /i^/3, where hv is the photon energy before the 
scattering. With typical CMB photon energy (hv ^ 2.3 x 
10 -4 eV), the range of electron's relativistic gamma fac- 
tor corresponding to the 10-40 keV inverse Compton hard 
X-ray becomes 3.3 x 10 7 < 7 2 < 1.3 x 10 8 . On the other 
hand, a typical synchrotron radio frequency (or, a criti- 
cal frequency) emitted by the electrons with the magnetic 
field strength B is (^ c /MHz) ~ 3.37 2 (£/G), where homo- 
geneous pitch angle distribution is assumed. Therefore, 
typical synchrotron radio frequency range emitted by the 
electrons in the above-mentioned energy range becomes 
1.1 x 10 8 (£/G) < (v c /MRz) < 4.3 x 10 8 (£/G). It is re- 
ported that radio flux of A2319 halo is 1.0 Jy at 610MHz 
and that its spectral index is 0.92 (408-610MHz) or 1.4 
(26-610MHz) (Harris & Miley 1978; Feretti et al. 1997). 
Thus, radio flux corresponding to the energy range of 10- 
40 keV inverse Compton hard X-ray becomes values as 
in table 7 for each spectral index, with a monochromatic 
approximation for a single electron's spectrum and on the 
assumption of a single power-law electron energy distri- 
bution. 

Flux of the inverse Compton scattering of CMB pho- 
tons and synchrotron radiation from the same electron 
population has the following relationship, 

Fic Ucmb Ucmb ^ 



mag 



B 2 /8w' 
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Table 7. Radio flux in the energy band corresponding to the 
10-40 keV inverse Compton component 



photon index radio flux (P S ync) 

(erg s _1 cm -2 ) 
L92 2.3 x 10- i5 (£/G)- u - U8 

2A 5.1 x lQ- 16 (ff/G)-°- 4 



Table 8. The lower limit of the magnetic field strength for 
each model. 



model B(/iG) 

lfcT + PL(1.92) >0.19 

l/cT + PL(2.4) >0.27 

2fcT + PL(1.92) >0.14 

2feT + PL(2.4) >0.25 



where Uqmb and J7 mag are the energy density of CMB 
photons and magnetic field, respectively. With Uqmb = 
5.2 x 10 -13 erg cm -3 at z = 0.0557 and the values in table 
6 and 7, the lower limit of the magnetic field strength 
for each model is obtained as in table 8. Although the 
detailed results slightly depend on the spectral modeling, 
basically the field strength tends to be more than ~ 0.2 
/iG. 

A2319 was observed by Beppo-SAX PDS in the hard 
X-ray band (Molendi et al. 1999). The obtained upper 
limit flux of the power-law component in the energy band 
of 13-50 keV is 2.3 xlO" 11 or 2.0 xlO" 11 erg s" 1 cm" 2 , 
depending on the adopted radio spectral index 0.92 (408- 
610 MHz) or 2.2 (610-1420 MHz), respectively. Different 
radio spectral indexes are used because the PDS results 
rely on higher energy band than PIN. When these results 
are compared with ours, however, it should be noted that 
the Crab's 20-80 keV flux derived from Beppo-SAX PDS 
is 21 % smaller than that derived from Suzaku. Thus, 
the above-mentioned Beppo-SAX results can be converted 
into Suzaku-equivalent fluxes of 2.9 xlO -11 or 2.5 xlO -11 
erg s _1 cm -2 for the radio spectral index 0.92 or 2.2, re- 
spectively. These results seem to be similar to our results 
of the lkT + PL(1.92) model at first glance. However, 
please remember that HXD PIN FOV is much narrower 
than Beppo-SAX PDS, which is preferable in order to 
avoid contamination sources, but large enough to cover 
the A2319 radio halo entirely. 

Recently, Ajello et al. (2009) studied the hard X-ray 
properties of 10 galaxy clusters including Abell 2319 us- 
ing Swift-BAT. The obtained upper limit of the power- 
law component in the energy band of 10-40 keV is 2.9 
xlO -12 or 1.7 xlO -12 erg s _1 cm -2 with a radio spec- 
tral index 0.92 through Swift-BAT alone or Swift-BAT + 
XMM-Newton analysis, respectively. Here, their original 
values at the 3a confidence level in the energy band of 50- 
100 keV are converted into those at the 90% confidence 
level in the 10-40 keV band. Considering that A2319 has 
relatively high temperature (~ 10 keV) and that Swift- 
BAT, with its long exposure, has better sensitivity in the 



higher energy band (~ 100 keV) than Suzaku PIN, it is 
not so surprising that their results give tighter constraint. 
Suzaku is superior in determining the thermal properties, 
thanks to its high sensitivity around 20-40 keV. We note 
that their power-law "upper-limit" is quite low, much less 
than those inferred from the sensitivity. Their sensitivity 
plot shows ~ 1 mCrab, or ~ 1.6 x 10~ n at 10-40 keV, as 
a 3-sigma sensitivity, or even worse above 100 keV. Thus, 
Swift-BAT is presenting ~ 2 or less xlO -12 erg s _1 cm -2 
upper-limit using the detector with a 90% confidence sen- 
sitivity of ^ 7 xlO -12 erg s _1 cm -2 . This is of course al- 
lowed in some cases, especially when the photon statistics 
favors negative flux. Anyway, our results, independently 
obtained, are consistent with theirs. We would also like 
to point out that the Swift alone analysis lacks soft X- 
ray spectral information important to determine thermal 
component precisely, and that generally speaking cross- 
calibration between detectors onboard different satellites 
(such as XMM and BAT) is not so easy. 

As we wrote before, ASCA results (Markevitch 1996) 
suggest there is lower temperature (~ 8 keV) ICM in the 
outer region without XIS but within PIN FOV, though 
this is not clear in Beppo-SAX results (Molendi et al. 
1999). The existence of this cold ICM could affect the 
estimation of the upper limit of the power-law compo- 
nent. To evaluate this, we make the spectra assuming 
that the outer region (r > 10') has 8 keV and perform 
joint spectrum analysis of XIS and PIN in a similar way 
in subsection 3.3. The resultant flux of the power-law 
component in 10-40 keV increases by 10.4 x 10 -12 and 
9.9 x 10 -12 erg s _1 cm -2 for the lkT + PL models with a 
photon index 1.92 and 2.4, respectively. These values are 
smaller than both the statistical errors and those induced 
by background uncertainty, as listed in table 6. This effect 
shall affect the Beppo-SAX PDS results as well. 

4-4- The Energy Budget of the Intracluster Space 

From the observed results, we can get information about 
the energy densities of the thermal ICM, non-thermal elec- 
trons, and magnetic field in the intracluster space of the 
A2319, which are basic and important parameters for the 
cluster structure and evolution. For simplicity, we as- 
sume that the radio halo region is a sphere with a ra- 
dius 0.66 Mpc. With /^-model parameters obtained by 
ROSAT (Trevese et al. 2000) and kT = 10 keV, energy 
density of the thermal ICM becomes /7th = 4 eV cm -3 . 
From the upper limit of an inverse Compton component 
of the lkT + PL(1.92) model in table 6, we obtain the 
energy density of the relativistic electrons corresponding 
to 10-40 keV hard X-ray (or, 5.7 x 10 3 < 7 < 1.1 x 10 4 ) 
U e <2x 10~ 2 eV cm" 3 . Therefore, U e /U th < 5 x 10~ 4 , 
though in this calculation we do not consider the contribu- 
tion from relatively low energy electrons, which could be 
dominant in the energy density of the non-thermal elec- 
tron populations. With the lower limit of the magnetic 
field for the lkT + PL (1.92) model in table 8, the energy 
density of the magnetic field becomes U mag > 1 x 10 -3 eV 
cm -3 , which means Z7 mag /Z7th > 3 x 10 -5 . 



Suzaku Observation of Abell 2319 



11 



4.5. A Particle Acceleration Scenario 

In this subsection, let us discuss a particle acceleration 
scenario of the A2319 implied by the XIS and HXD results 
mentioned above. The XIS results tell us two important 
points as follows. First, the central cool region associated 
with the cold front have higher metal abundance, which 
suggests that this might be an old remnant of a subclus- 
ter's cool core that infalls in the past. Second, ICM ve- 
locity difference in the observed region is probably sub- 
sonic, which means that no definitive shocks are expected 
there. This is also supported by the fact that there is 
not a significant very hot component from the XIS+HXD 
analysis. On the other hand, existence of subsonic tur- 
bulence (Av ~ 500 km s _1 ) can be allowed. Therefore, it 
is more likely that non-thermal electrons relevant to the 
radio halo of the A2319 are accelerated by the intraclus- 
ter turbulence rather than the shocks. The turbulence 
is probably excited by the past merger activity, which is 
also responsible for the cold front and central cool core 
with higher abundance. This picture is also consistent 
with the fact that numerical simulations show that the 
turbulence motion is developed in the late phase of merg- 
ers (Takizawa 2005). The discussion on the energy budget 
in the previous subsection shows that the magnetic field 
energy density is likely much less than the hydrodynami- 
cal turbulence, which means that the magnetic turbulence 
can be easily excited (Roettiger et al. 1999; Asai et al. 
2007; Takizawa 2008). Although a collision between the 
subgroup A2319A and A2319B is likely supersonic, it does 
not seem to form sufficiently strong shocks to excite the 
radio halo activity. It is probable that they are not so 
close to each other that sufficiently strong shocks occur. 

4-6. Future Prospects 

We obtained only an upper limit of the ICM line-of-sight 
velocity difference. To resolve the intracluster turbulence 
(Av ~ 500 km s _1 ), a different type of X-ray spectrome- 
ter than CCD is highly desired. The ASTRO-H satellite 
(Takahashi et al. 2006), which is planned to be launched 
around 2013, will enable us to measure the line-of-sight ve- 
locity directly with the X-ray microcalorimeters and pro- 
vide us with useful information on the dynamical status 
of the ICM. ASTRO-H is also planed to have hard X-ray 
mirrors and imager, which enable us to get hard X-ray 
images at very high sensitivity. This gives us very im- 
portant clues to understand where and how the particle 
acceleration processes occur in the intracluster space. 

5. Conclusions 

We observed the central region of the Abell 2319 with 
Suzaku. From the XIS analysis, we found that the cen- 
tral cool region in the north-west of the cold front has 
higher metal abundance than the surrounding region for 
the first time. We measure the line-of-sight velocities of 
the ICM with XIS and find that the velocity difference 
is less than ~ 2000 km s _1 , which means that the ICM 
motion is probably subsonic. In addition, we show that 



the velocities in the observed region are consistent with 
those of entire A2319 and A2319A subgroup for the first 
time. No signs of ICM motion related to A2319B subgroup 
are found. If the ICM associated with the A2319B sub- 
group has the similar velocity of A2319B member galax- 
ies, a supersonic collision is expected in the ICM. From 
the XIS+HXD wide band spectral analysis, we search for 
very hot and non- thermal inverse Compton components. 
There can be a little amount of very hot (~ 15 keV) com- 
ponent, but it does not have to exist. No significant in- 
verse Compton component is detected, and we derive its 
upper limit for several spectral models. Typically, it be- 
comes < 3 x 10 -11 erg s _1 cm -2 , which is consistent with 
the recent Swift-BAT results. Comparing the synchrotron 
radio flux and upper limit of the inverse Compton one, 
we obtain the lower limit of the magnetic field strength 
B > 0.2/iG, though the results slightly depend upon the 
spectral modeling. Taking into account of the lack of the 
supersonic ICM motion and a very hot component, we 
conclude that the relativistic electrons responsible for the 
radio halo of the Abell 2319 are more likely accelerated 
by the turbulence rather than shocks. 
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